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In the context of open inflation, we calculate the probability distri- 
bution for the density parameter £1. A large class of two field models of 
open inflation do not lead to infinite open universes, but to an ensemble 
of inflating islands of finite size, or "quasi-open" universes, where the den- 
sity parameter takes a range of values. Assuming we are typical observers, 
the models make definite predictions for the value f2 we are most likely 
to observe. When compared with observations, these predictions can be 
used to constrain the parameters of the models. We also argue that ob- 
severs should not be surprised to find themselves living at the time when 
curvature is about to dominate. 



I. INTRODUCTION 

Anthropic considerations have often been used in order to justify the "naturalness" of 
the values taken by certain constants of Nature ffl. In these approaches, it is assumed that 
the "constants" are really random variables whose range and "a priori" probabilities are 
determined by the laws of Physics. Knowledge of these "a priori" probabilities is certainly 
useful, but not sufficient to determine the probability for an observer to measure given 
values of the constants. For instance, some values which are in the "a priori" allowed 
range may be incompatible with the very existence of observers, and in this case they will 
never be measured. The relevant question is then how to assign a weight to this selection 
effect. 

A natural framework where these ideas can be applied is inflation. There, the false- 
vacuum energy of the scalar field which drives the inflationary phase can thermalize in 
different local minima of its potential, and each local minimum may have a different set 
of values for the constants of Nature. Also, there may be different routes from false 
vacuum to a given minimum. In this case all thermalized regions will have the same low 
energy Physics constants, but each route will yield a hot universe with different large 
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scale properties. Here, we shall be concerned with this possibility, where the fundamental 
constants (such as the gauge couplings or the cosmological constant) are fixed, but other 
cosmological parameters such as the density parameter or the amplitude of cosmological 
perturbations are random variables whose distribution is dynamically determined. 

In this context, the most reasonable -and predictive- version of the anthropic prin- 
ciple seems to be the principle of mediocrity according to which we are typical 
observers who shall observe what the vast majority of observers would. Thus, the mea- 
sure of probability for a given set of constants is simply proportional to the total number 
of civilizations emerging with those values of the constants. In this paper we shall use 
this principle in order to calculate the probability distribution for the density parameter 

n. 

Standard inflationary models predict Q = 1 with "certainty". What this means is 
that these models can explain the observed homogeneity and isotropy of the universe only 
if the universe is flat. However, a class of "open inflation" models which lead to fi < 1 
have received some attention in recent years |7]-|9|]. In these models, inflation proceeds 
in two steps. One starts with a scalar field a trapped in a metastable minimum of its 
potential V(a). The false vacuum energy drives an initial period of exponential expansion, 
and decays through quantum nucleation of highly symmetric bubbles of true vacuum. 
The interior of these bubbles has the geometry of an open Friedmann-Robertson- Walker 
universe. This accounts for the observed homogeneity and isotropy of the universe. In 
order to solve the flatness problem a second stage of slow roll inflation inside the bubble 
is necessary. 

In models with a single scalar field a, all bubbles have the same value of Q which is 
determined by the number of e-foldings in the second period of inflation. The potential 
V(a) in such models is assumed to have a rather special form, with a sharp barrier next 
to a flat slow-roll region, which requires a substantial amount of fine-tuning. Additional 
tunning is needed to arrange the desired value of Q. A more natural class of models 
includes two fields, a and 0, with a doing the tunneling and <f> the slow roll ||. The 
simplest example is 

V(a,<P) = V t (a) + ^ 2 <P\ (1) 

where Vo(cr) has a metastable false vacuum at a — 0. After a tunnels to its true minimum 
a = v, the field would drive a second period of slow roll inflation inside the bubble. 
Depending on the value of <fi at the time of nucleation, the number of e-foldings of the 
second stage of inflation would be different. 

Initially, it was believed |§ that models such as (fj) would yield an ensemble of 
infinite open universes, one inside each nucleated bubble, and each one with a different 



value of the density parameter. However, it has been recently realized Jl0| that this 
picture is oversimplified. The two field models which allow for variable Q do not actually 
lead to infinite open universes, but to an ensemble of inflating islands of finite size inside 
of each bubble. These islands are called quasi-open universes. Within each island, the 
number of e-foldings of inflation decreases as we move from the center to the edges. Also, 
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each island is characterized by a different number of e-foldings in its central region. As 
a result, even within the same bubble, different observers will measure a range of values 
of the density parameter. The picture of the large scale structure of the universe in these 
models is rather simple, because all bubbles have the same statistical properties. We shall 
see that the quasiopen nature of inflation is of crucial importance for the calculation of 
the probability distribution for the density parameter. 

In models of quasiopen inflation, such as ([!]), Q takes different values in different 
parts of the universe, while the other constants of Nature and cosmological parameters 
remain fixed. More general models can be constructed where other parameters can change 
as well, and in Section VII we give an example of a model with a variable amplitude of 
density fluctuations. However, our main focus in this paper is on the models in which 
only Q is allowed to vary. 

In order to apply the principle of mediocrity to our models, we will have to compare 
the number of civilizations in parts of the universe with different values of Q. Of course, we 
cannot calculate the number of civilizations. However, since the value of Q does not affect 
the physical precesses involved in the evolution of life, this number must be proportional 
to the number of habitable stars or, as a rough approximation, to the number of galaxies. 
Hence, we shall set the probability for us to observe a certain value of Q to be proportional 
to the number of galaxies formed in parts of the universe where Q takes the specified value. 

The principle of mediocrity was applied to calculate the probability distribution for 
Q in an earlier paper Jy]] , which assumed the old picture of homogeneous open universes 
inside bubbles. A serious difficulty encountered in that calculation was that open universes 
inside the bubbles have infinite volume and contain an infinite number of galaxies. Thus, 
to find the relative probability for different values of Q, one had to compare infinities, 
which is an inherently ambiguous task. This problem was addressed in [|ll]] by introduc- 
ing a cutoff and counting only galaxies formed prior to the cutoff. Although the cutoff 
procedure employed in |Tl| has some nice properties, it is not unique, and the resulting 



probability distribution is sensitive to the choice of cutoff [ 12| . This cutoff dependence 



which also appears in other models of eternal inflation [12,13 1 , has lead some authors to 



doubt that a meaningful definition of probabilities in such models is even in principle 
possible 00. 

However, this pessimistic conclusion may have been premature. According to the 
quasiopen picture, fl takes all its possible values within each bubble. Since all bubbles 
are statistically equivalent, it is sufficient to consider a single bubble. Moreover, we can 
restrict ourselves to a finite (but very large) comoving volume within that bubble, provided 
that its size is much greater than the characteristic scale of variation of Q. Thus, we no 
longer need to compare infinities, and the problem becomes well defined. 

The possibility of unambiguous calculation of probabilities in the quasiopen model 
was our main motivation for revising the analysis of Ref. ||11|| . Also, we shall give a more 
careful treatment of the astrophysical aspects of the problem which were discussed rather 



sketchily in [11 



The paper is organized as follows. In Section II we review the main features of quasi- 
open inflation. In Section III we introduce the probability distribution for Q. A basic 
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ingredient in this distribution will be the anthropic factor which gives the number 

of civilizations that develop per unit thermalized volume in a region characterized by a 
certain value of Q. In Section IV we evaluate and calculate the probability distribu- 
tion for Q in the model (|1]). In Section V we extend our results to more general models 
with arbitrary slow roll potentials for the field <fi. In Section VI we discuss observational 
constraints on quasiopen models due to CMB anisotropies and how these constraints re- 
strict the class of models that give a probability distribution peaked at a non-trivial value 
of Q. In Section VII we comment on the "cosmic age coincidence" , that is, on whether it 
would be surprising to find ourselves living at the time when the curvature of the universe 
starts dominating. In Section VIII we summarize our conclusions. Some side issues and 
technical details are discussed in the appendices. 



II. QUASI-OPEN INFLATION 

In this Section we shall review the main features of quasi-open models which will 
be relevant to our discussion. To begin with, we shall consider a model of the form ([[]). 
In Section V we shall consider more general slow-roll potentials. 

As mentioned in the introduction, the interior of a bubble is isometric to an open 
Friedmann-Robertson- Walker universe, with line element 

ds 2 = -dt 2 + a 2 {t) [dr 2 + sinh 2 r{dd 2 + sin 2 6dtp 2 )}. (2) 

The scale factor a obeys the Friedmann equation 

o fa\ 2 8ttG 1 , . 

\aj 6 a z 

At sufficiently early times (t — > 0), the curvature term in the r.h.s. dominates over the 
energy density p of the scalar fields, and the scale factor behaves as a « t. 

For the second period of slow roll inflation inside the bubble, the energy density of 
the scalar fields must be dominated by the potential term 

y(a,0)» ( r 2 ,0 2 . (4) 

Inside the bubble, the field a quickly settles down to its v.e.v. a = v with V t (v) = 0. 
This happens on a time scale of order t ~ M _1 , where M is the typical mass scale of V t . 
After that, becomes a free field with constant mass 

2 9 2 

m — —v , 
2 

and the condition for inflation becomes 0~M P , where M 2 = G~ x . 

An important feature of quasi-open models is the existence of the so-called super- 
curvature modes for the slow roll field 0. These are modes which are not normalizable on 
the infinite t = const, hyperboloids inside the bubble, but which nevertheless have to be 
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included in the field expansion. The reason is that they are normalizable on the Cauchy 
surface where equal time commutation relations are imposed. Supercurvature modes are 
characterized by their eigenvalue of the Laplacian on the 3-hyperboloid. For the model 
(P, this eigenvalue is given by ]TU 



7 = « 1. (5) 

Here, 

Hi = S -fv o (0) 

is the Hubble rate during the first stage of false vacuum dominated inflation and Rq is 
the size of the bubble at the time of nucleation, which can be given in terms of the model 
parameters Typically, m 2 <C Hp < Rq 2 , from where the condition 7 <C 1 follows. 



Around the time to when a settles down to its v.e.v., the field <p will be in a homo- 
geneous and isotropic quantum state with mean squared amplitude given by |L0| 

'"-<*>"'(§)? (6) 

The presence of the factor 7 _1 reflects the fact that / is dominated by the contribution 
of supercurvature modes. Introducing 7 from (|), we find 

7T TuRq 

Up to numerical factors, this is basically the finite temperature dispersion of a field of 
mass m at the Rindler temperature given by T = (27i\Ro) _1 - The correlation length of (f) 
is given by [[jlj r ~ 7 -1 . This means that at the time t ~ i we can divide the space into 
regions of co-moving size r ~ 7 _1 ^> 1 where the field is coherent. Notice that the size of 
these regions is much larger than the curvature scale r — 1. 

The parameters of the model can be chosen in such a way that / is close to the 
Planck scale, and in that case the slow roll field easily reaches inflating values <fi ~ Mp. 
For instance, if the potential V t is such that the bubble walls are thick, then R Q ~ M~ l . 
Taking M ~ 10 16 GeV and m ~ 10 13 GeV, we find / ~ M p . In this case, inflating regions 
of co-moving size r ~ 7 _1 where the field is large and positive will be next to inflating 
regions where the field is large and negative. These two inflating regions will be separated 
by regions where the field is small and the universe does not inflate. 

The parameters can also be such that / <C M p , and in that case most of the regions 
will not attain an inflating value of 0. Inflation will only happen in those regions where, 
as a result of a statistical fluctuation, the field happens to be far above its r.m.s. value. 
Since the volume of the hyperboloid is infinite, there will be a small but finite density of 
these inflating islands inside of each bubble. Those rare "high peaks" will have spherical 
symmetry. If we take the inflating patch to be centered at r = 0, the radial profile of the 
field is given by 
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,. . , sinh[(l — 7) 1 / 2 rl 
0(t o ,r)^0 o ' 



'1 — 7) 1 / 2 sinhr ' 

Where 0o = <K^o ; 0) is a constant. The probability distribution for O is given by 



P(<f>o) oc exp 



2p 



(9) 



The variation of 0o within the bubble results in a position-dependent number of inflation- 
ary e-foldings, and thus in a variable density parameter Q. Note that all other cosmological 
parameters, such as the amplitude of the density fluctuations, remain fixed throughout 
the bubble (and are the same for all bubbles). The probability @ is one of the basic 
ingredients from which the most probable value of Q is calculated. 

It should be mentioned that the size of an inflating region can be much larger 
than the size of the actual "populated" region within it, r p , where matter will cluster 
efficiently into galaxies. The size of the populated region is calculated in Appendix A. 
This size should be larger than the present horizon, since otherwise we would observe large 
anisotropics in the galaxy distribution. For Q not too close to 1, the horizon distance is 
comparable to the curvature scale, r = 1, and we have to require that r p > 1. The 
corresponding constraint on 7 is obtained in Appendix A. 

Equation @ can be understood from a different perspective, by using the Euclidean 
approach to the calculation of the nucleation rate. The strategy is to study how this rate 
is affected by the local value of at the place where the bubble nucleates. This is simple 
because we only need this Euclidean action to quadratic order. Taking = 0, we denote 
by o"o( r ) the 0(4) symmetric instanton responsible for vacuum decay. Here, r = it is 
the Euclidean time, the "radial" coordinate on which the instanton depends. Expanding 
the Euclidean action Se to second order in perturbations of a and of 0, the perturbations 
in cr and will decouple to quadratic order. Taking = O = const, the change in the 
Euclidean action will simply be 

AS E = J (g/2)al(T)<Pld A x (10) 

We can approximate the integral by taking ga\ = m 2 inside the volume of the bubble and 
(To = outside. Then we have 

~2 u 2 

From the formula P ~ exp(— Se), we essentially recover Equation (|9|). Even though 
we have used the thin wall approximation, we should stress that the coincidence of this 
"adiabatic" result (where the field <ft is taken as constant) with the field theoretic one 
(where <p is quantized in the bubble background and its r.m.s. is evaluated right after 
bubble nucleation) is also valid for thick walls JT0[ . 

As emphasized in [|l(]], the adiabatic approach to the calculation of the distribution 
of 0o should be interpreted with caution. It does not mean that the surface t = to inside 
the bubble will have a constant value of the field 0. It only gives the probability that a 
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ASe = -m 2 R 4 ^ 
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bubble will nucleate with the value of 4> = 4>q near r = 0. We know, however, that the 
quantum state of a nucleating bubble is homogeneous, and therefore in the ensemble of 
bubbles there is nothing special about the point r = 0. Therefore, this also gives the 
probability distribution for around any point inside the bubble. 



III. PROBABILITY DISTRIBUTION FOR O 



In this Section we shall follow some of the steps used in Ref. |TT| for the calculation 



of the probability distribution for f2, although the present case will actually be simpler. 



In the case of Ref. ||11|| , one had to deal with an infinite number of bubbles, each one con- 
taining an infinite open universe with a different density parameter. Since the probability 
for a given set of parameters is roughly proportional to the total volume that ends up 
having those values of the parameters, one had to face the difficulty of comparing infinite 
volumes in an eternal inflationary universe . 

In our case, all bubbles are statistically equivalent. All of them are described by a 
homogeneous and isotropic quantum state, with (4> 2 ) given by (|6|). Hence, in order to find 
the probability distribution for O it is sufficient to look at the interior of a single bubble. 
Also, since the quantum state is homogeneous, we only need to consider the evolution of 
a patch of finite co-moving size around an arbitrary point on the t ~ M _1 hyperboloid. 
The patch should be sufficiently large that it contains regions with all possible values 
of (ft, distributed according to (|9]). Since inflation inside the bubble is not eternal, the 
number of civilizations resulting from this co-moving patch is finite and there is no need 
for regularization. 

As mentioned in the previous section, at early times the scale factor behaves as 
a ~ t. By the time t ~ H~ 1 (<po), where 

H(<j) ) = (47rG/3) 1/2 m0 o , 

the energy density in the scalar field <j) starts dominating over the curvature term. If the 
condition for slow roll inflation 

O > <j> th = -^L (11) 

V47T 



is satisfied, then using equation (^) the scale factor will subsequently evolve as 

a(t) » H^e»®, 

where H = H(cf) ) and 

N(t)- 

Using the slow-roll equation of motion for (f) 

2 

• — m 

= w 



m H{(j>) 
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we have N(t) ~ 2tiG [4>q — (fi 2 (t)}. Since 4>o is actually a slowly varying function of position, 
the scale factor is a local one, and should be understood as a(t,x l ). Notice that the co- 
moving scale over which a changes is comparable to 7 -1 and hence it is much larger than 
the curvature scale, so it is meaningful to use the Friedmann equation (|]). 
The number of e-foldings of inflation depends on the local value of 0o 



(hhVPo) 



#-i e "tfc(A0 ~ tf-V^o-^, (12) 



where (p t h is defined in (|TT|). It will be convenient, as a first step, to find the probability 
distribution for a random "civilization" to live in a region which had a value of the slow 
roll field equal to 0o & t the beginning of inflation. This is given by 

dV((f>o) = P(0oK\(0oM0o)d0o. (13) 

Here 

P(0o)oce2/ 2 (14) 

is the probability that a given point on the t ~ M _1 hyperboloid will have the value 0o 
right after nucleation. Because the number of civilizations is proportional to the volume, 
we have inserted the total expansion factor during inflation a\ h {4>o)- Finally, u(<po) is the 
"human factor", which represents the number density of civilizations that will develop 
per unit thermalized volume function of 0o- 
Factoring out the dependence on u, 

dV = is((f>o)dV, (15) 

the leading exponential behaviour of dV is 

1 



exp 6nG — 



2/2 



where we have used ( 12|) and (|13|). The behaviour of V depends on whether / is large or 



small compared with M p . Defining 

^^f' (16) 

it is clear that for fi < 1/2 large values of 0o are favoured due to the gain in volume factor, 
and we may expect the universe to be very flat. For fj, > 1/2, the volume factor alone is 
not sufficient to compensate for the exponential suppression of high peaks. We shall see 
that the human factor may play an important role in this case . 

It is convenient to express the above distribution in terms of the density parameter. 
Following 0,11 1, we have 



[HiMa^o)] 2 = 1 + B-^ « B^. (17) 



S 



Here 

T 2 

B ps — — — , 

TeqTcMB 

T eq is the temperature at equal matter and radiation density, T t h is the thermalization 
temperature and Tqmb is the temperature of the cosmic microwave background, measured 
at the same time as Q. Typically, B is exponentially large, with (In B) ~ 10 2 . ^From (|T3| ) 
and (|T7]) we find 

/ s ? f d\na f h\ 1 dtt . . 

dV ^ = p ^<{-l^) mTny < 18 > 

Using (0), and disregarding the logarithmic dependence on Q, we find 

dP (SI) oc n 1 ' 2 - 3 ^! - fi) 3 ^ 5 / 2 c/f2, (19) 



where fx is given by Eq. (|16|). For /i > 5/6 the probability distribution is peaked at f2 = 0, 
for /i < 1/6 it is peaked at Q = 1, and for the intermediate range 1/6 < /x < 5/6 it has 
two peaks, one at Q = and one at O = 1. However, it is easily seen that for \i > 1/2 
the highest peak will be at Q = whereas for [i < 1/2 it will be at O = 1.Q Note that all 
dependence on the particle physics model in Eq.(|T9"D has been compressed into a single 
parameter li. 

Eq. (TOT) is the same expression that was found in Ref. JTT1] by considering an en- 



semble of bubbles with different values of fl and using the prescription introduced in 
fl"3H for the regularization of infinite volumes. We regard the agreement between the two 
approaches as a validation of this regularization prescription (in models where regulariza- 
tion is needed). Alternative regularizations proposed in |12| give different results and are 
therefore disfavoured. 

Let us now include the human factor u(Q). As mentioned above, this will play a role 
for /i > 1/2, when the expansion alone is not sufficient to compensate for the exponential 
suppression in O due to tunneling. Since the probability distribution V tends to peak 
near the extremes, it is convenient to work with a logarithmic variable which gives equal 
measure to each decade in the vicinity of Q — or Q — 1. One such variable is lnx, where 

i-n 

x - — • (20) 
Hence, we shall be interested in the probability density 

W(Q) = « v(x) x^' 1/2 \ (21) 
dmx 



Strictly speaking, the peak would not be exactly at 0, = 1 because the Gaussian distribution 
dl4| ) is only an approximation which ignores the backreaction of the slow roll field on the bubble 
background. We shall return to this issue in Section V. 
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The peak of this distribution will give the most probable value of £1 

It should be noted that, since the density parameter changes with time, both v(Q) 
and iy(£2) are in principle time dependent. However, this time dependence is somewhat 
trivial, entering (^1|) through the parameter Tqmb, the temperature at which the density 
parameter is equal to Q. What we are actually interested in is the probability distribu- 
tion for different types of thermalized regions, which is intrinsically time independent. We 
could, for instance, set T C mb equal to the temperature at recombination, and then the 
probability distribution would be expressed in terms of f2 rec , which completely character- 
izes the history of a given thermalized region. Noting that Friedmann's equation can be 
rewritten as x~ x = (87rG/3)pa 2 , in the matter era we have 

x oc a(t) oc Tcl IB , (22) 

where a indicates the scale factor. Hence, in practice, we can use a "gauge invariant" 
approach: we shall write W(fl) as a function of the product x Tcmb, which is time 
independent in the matter era. 



IV. THE ANTHROPIC FACTOR 

In previous work [p|- |i9| . |iT| . |2~0f , was taken to be proportional to the fraction of 
clustered matter f c on a relevant mass scale M g . This scale can be chosen as the typical 
mass of an galaxy, M g ~ 1O 12 M |I7| , |20[] , given that most of the observed luminous 



matter is in this form. Also, galaxies much smaller than 1O 12 M may not be suitable for 
life, because their gravitational potential would not be able to hold the heavy elements 
produced in supernovae explosions. Matter will only cluster when the density contrast 
S(M g ) extrapolated from linear perturbation theory exceeds a certain threshold 5 C . Hence 



f c can be estimated as |2T|j2~ 



f c (M g ,t)=erfc(-^-). (23) 



y/2, 



a . 



Here erfc is the complementary error function and cr(M g , t) is the dispersion in the density 
contrast, also evolved according to linear theory Q |fj3| 

a(M 9 ,t) = %W /( x), (24) 



2 This expression for the growth of perturbations is different from the one used in [11|. There, 
the growth factor from the time of equilibrium of matter and radiation was considered, and a 
spurious factor of O was included, which was actually due to the uncertainty in the value of the 
redshift at the time of equilibrium. This factor should actually not be present in the probability 
distribution for f2, since the time of equilibrium is the same in all thermalized regions. 
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where x is given by (pPf) and 



f(x) = 1 + - + 3(1 + ! )V2 ln[(l + xf* - x 1 / 2 }. (25) 



x x 3 / 2 

The subindex rec denotes quantities evaluated at the time of recombination. In an open 
universe, perturbations stop growing after the universe becomes curvature dominated. 
Since we are interested in the total fraction of clustered matter in the entire history of a 
given region, we should use in (|2"3"| ) the asymptotic value of a at large times (x — > oo), 
which approaches a constant. 

In a flat universe, the critical density contrast takes the value 5 C ~ 1.7. However, it 
is known that S c should be slightly ^-dependent [p^j . The variation is rather small, and 



5 C changes by no more than 5% as VL varies from 0.1 to 1. Here we adopt the value of 5 C 
estimated in the spherical collapse model as p5|j26|| 



6 c (x) = \f{x)g{x\ (26) 

where 

( 7T \ 2/3 

9{x) = 1 + Uv2(i + x y/2 _ sinh -i x i/2 J • ( 2? ) 

For x — > 0, we have 5 C = (3/5)(37r/2) 2//3 w 1.69, as in the case of a flat universe, and for 
x — >• oo we have <5 C = 3/2. 

Substituting (|2~4j) and (^) in Eq. (^) and taking the limit x — ► oo we obtain 

^ = erfc iy ^J^ ^ = erfc(y). (28) 

The distribution ( P5| ) is given as a function of the density parameter at the time of re- 
combination. As mentioned at the end of the last section, in order to compare predictions 
with observations, it is convenient to express the distribution as a function of x at any 
temperature T C mb- Using (|22| ) we have 

_ 3 T C mb 
y = kx = — -= — x. (29) 

0\ Z(T rec J- rec 

In order to evaluate the coefficient k, we need to know a rec . It is clear that a rec has 
nearly the same value in all regions where curvature dominates only well after the time of 
recombination. In principle this value is given in terms of the parameters of our theory 
of initial conditions. 
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FIG. 1. The coefficient k which relates the variable y to the density parameter 
y = kx = k(1 — f2)/0, depends on a rec , the value of the density contrast at the time of re- 
combination. Our ability to infer a rec from CMB observations is limited by the fact that f^o 
in our observable "subuniverse" is not known very precisely. In the figure we plot the inferred 
value of k for various assumed values of Qq. The value of a rec depends moreover on the scale 
Rgai corresponding to objects of galactic mass. The curve is plotted for two different values of 
this scale (see Appendix B). The parameter k depends on the temperature at which we observe 
fi. Here we have taken Tqmb = 2.7K . 

In practice, we can adjust the parameters of the theory to fit CMB observations. Our 
ability to infer a rec from CMB observations is, however, limited by the fact that it depends 
on the values of Qo and h in our visible universe, which are not very well determined. As 



noted in [20], this limitation also arises in attempts to find the probability distribution for 
the cosmological constant. Therefore, until determinations of o rec become more precise, 
the best one can do is to assume certain values of Qq and h and check whether the assumed 
values fall within the range favoured by the resulting probability distribution for Q. The 
value of k for Tqmb = 2.7 K is estimated in Appendix B and plotted in Fig. 1 as a 
function of Qq. For each value of Qq, h has been chosen so that the "shape parameter" 
T « Q h « .25 (see Appendix B). Also, there is some uncertainty in the relevant co- 
moving scale Rg a i corresponding to M ga i [p0[| . I n the figure we consider two possibilities, 
Rgai = 1 h~ x Mpc and R ga i = 2 h^Mpc. For fl in the range 0.1 < Q < 0.7 we find that 

K ~ 0.1. 
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0.01 
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FIG. 2. The probability distribution pl| ) as a function of y, for various values of /j,. Also 
represented is the fraction of clustered matter v{y) as a function of y. 

The fraction of clustered matter v is shown by a solid curve in Fig. 2 as a function 
of y. In our universe, the density of matter presently clustered in giant galaxies satisfies 
Qgai > 0.05 f27j, which implies u(y) > 0.05S1 1 . The asymptotic value u(y) should be 



even larger. Solving for y, we obtain the observational constraint 

y < 0.9 . (30) 

The distribution 

W(Q) = oc erfc(y)i/ 3{ ^ 1/2) (31) 

amy 

gives the probability that a randomly selected civilization is located in a region which 
had a specified value of Q at a given temperature T C mb- It is represented in Fig. 2, as a 
function of In?/, for different values of the parameter //. The peak value y pea k, found from 
dWf dy = is plotted in Fig. 3 as a function of /i (curve a) . 
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FIG. 3. Peak of the probability distribution ( pip (curve a). The approximate value of y pea k 
given by (^) is represented by the curve b. Curve c represents the possible effect of helium line 
cooling failure, as discussed in Appendix C. 

For 2/^1, the error function can be approximated by 

erfc(y) » -Le^, (32) 
and the peak value can be expressed analytically, 

This curve is also shown in Fig. 3 (curve 6). Eq. (|33| ) can be rewritten as 

J peak \2 4/ 

which gives the peak value for the density parameter at the temperature T CMB . To 
estimate the width of the distribution (|3T| ) we expand lnjy to quadratic order in Amy 
around y peak , 

W « W peafc exp [-(3/i - 5/2)(Alny) 2 . 

Hence, the root mean squared dispersion in Q around its peak value will be given by (for 
fi Z 3/2) 
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Aln (!^5)~(6 M - 5)- 1/2 , (35) 

while the dispersion in y is independent of fj,, Ay ~ 1/2. 

From Fig. 2, we see that as /J is increased, the probability distribution is sharper 
and displaced towards larger values of y, in agreement with ( |3"3"| ) and fl3"5]). For /x = 1, 
the distribution has a substantial overlap with the region where (|30"D is satisfied and the 
fraction of clustered matter is compatible with observations. For /i = 5/2 this overlap is 
smaller, but still non- negligible. However, for fi = 11/2, the probability density at the 
point y ~ 0.9 is more than two orders of magnitude smaller than at its peak value. Particle 
physics models which give such high values of fi are therefore disfavored by observations. 
As we discussed in Section III, the probability distribution is peaked at £1 = 1 for fi < 0.5. 
Hence the range of \i that is of interest to us in these paper is 

0.5 < n < 3 (36) 

This corresponds to 

< y P eak £ 2. (37) 

It should be noted that the peak value for the fraction of clustered matter v(y pe ak), depends 
only on /i, and not on the primordial spectrum of density fluctuations. 

So far we assumed that all galactic- size objects collapsing at any time will form 
luminous galaxies. However, this is not necessarily the case. Galaxies forming at later 
times have lower density and shallower potential wells. They are vulnerable to losing all 
their gas due to supernova explosions ||28|| . Moreover, a collapsing cloud will fragment into 



stars only if the cooling timescale of the cloud t coo i is smaller than the collapse timescale 
T grav, otherwise the cloud would stabilize into a pressure supported configuration [ffij| , p8|j . 
The cooling rate of such pressure supported clouds is exceedingly low, and it is possible 
that star formation in the relevant mass range will be suppressed in these clouds even 
when they eventually cool. Hence, it is conceivable that galaxies that fail to cool during 
the initial collapse give a negligible contribution to u(Q) ||28|| . The possible effect of 
cooling failure and related phenomena on the probability distribution for Q is discussed 
in Appendix C, where we show that the effect is to shift the peak of the distribution 
towards larger values of fl These effects may be significant, but not dramatic, and ( |3~3| ) 
remains valid by order of magnitude. As an illustration, curve c of Fig. 3 shows the peak 
of the modified distribution when matter which clusters after the time when helium line 
cooling becomes inefficient is excluded from the anthropic factor 



V. MORE GENERAL MODELS 

In this section, we shall generalize our results to models where the slow roll potential 
is not necessarily quadratic in </>. In this case, the factor P(<po) i n Eq. (|T3|) can be estimated 
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in the adiabatic approximation, where the field </> is treated as a constant during tunneling, 
as described at the end of Section II. In this approximation we have 



P(0o) oc e 



where Se{4>o) is the action of the instanton for bubble nucleation, with the slow roll field 
frozen to the value 0o- 

^From (|I5D, we have 

W = —— oc erfc(y) 2 /- 3 / 2 P(0o) J~\ (38) 
amy 

where, as before 

^ _ 5 C 2 Tqmb 

\p2(J rec 5 Tree ^ 

and we have used ( |TTD to express the scale factor as a function of y. The Jacobian J is 
given by 



J 



dlny 



where we have used = H(<f)Q)e N , and the relation between the hubble rate H((p ) and 
the slow roll potential V in true vacuum 

Here <jt is the value of the tunneling field in true vacuum. We have also used the slow 
roll expression for the number of e-foldings 



iV(0 o ) = 8ttG 



<t>o V 



Here, as in (PU[), V stands for the derivative of V with respect to the slow roll field. 
Introducing the slow roll parameter 

^=if-Mv) 2<<1 ' <4o) 

we have 

j«(1^) 1/2 . 

In many models, the parameter e hardly changes in the relevant range of </>o, and hence 
we shall treat it as a small constant parameter. 

Extremizing (|38|) , we find that the peak value of 0o is given by the condition 



1 1 dlnerfc(y) 

H{(P0)\peak 



2 3 dlny 
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(41) 

peak 



where 



Eq. (f41"|) is the same condition we found in Section IV, and which is plotted in Fig. 3, 
except that now fi is a function of <po, and hence of y. 

Before we proceed, let us go back to the case discussed in Section III of a free slow 
roll field . Strictly speaking, the expression (|i~4f) for P(4>q) is just an approximation which 
is valid only for sufficiently low <fro, when the backreaction of <f> on the bubble background 
can be neglected. Now we can take this effect into account. For definiteness, let us 
consider the case where m ~ 10 13 GeV, and where the tunneling potential Vt(a) is such 
that false and true vacuum are strongly non-degenerate when = 0. In this case, the 
radius of the bubble is Rq ~ M _1 (thick wall bubble), where M ~ 10 w GeV is a typical 
mass scale in the tunneling potential. Let us denote by 4>deg the value for which the energy 
density corresponding to the slow roll potential is equal to the false vacuum energy in the 
unbroken phase Vp = V t (a = 0), 

2 m <l>deg = V F ~ — , 

where A is a self-coupling of the tunneling potential. For O <C (fideg, the value of /i is 
almost independent of O (this is the situation considered in Section III) 



/i(0 O < 4>deg) ~ A*0 



7r m 2 
48GM 4 ' 



The masses M and m can be easily adjusted so that jiQ <C 1/2. However, for O ~ 4>deg, 
the Euclidean action Se(0o) increases very steeply with (f> , and so does In this case, 
the condition (f4T|) will be satisfied for O ~ 4>deg ~ ^(A/Uo)" 1 / 2 , where M p = G^ 1 ^ 2 
is the Planck mass. The corresponding number of e-foldings of inflation is given by 
iV(0 o )^2vrG0 2 e ^(7r 2 /12)(A /U o)- 1 . ' 

Therefore, for /i 1/2, and with a suitable choice of A, the peak in the distribution 
may be adjusted to correspond to iV ps (1/2) In 5 w 60, where 5 is defined in (|T^). This 
is compatible with an open universe. However, this case is somewhat trivial, in the sense 
that the universe can be open only if the maximum allowed value of the slow roll field 
after tunneling, 0o = <fideg, does not drive a long enough period of inflation to make it flat. 

Turning to the general case, a more interesting situation arises when /i(0o) > 1/2 
throughout the range of <po (see Section III). In this case the product of tunneling and 
volume factors would peak at 4>o = 0, where the resulting universe would be almost empty, 
and the anthropic factor z/(f2) is crucial in determining the probability distribution for Q. 



3 Indeed, as approaches <fideg, the thin wall approximation starts to apply. Then, from Eqs 
( |43| ) and (f5|) below, we find that the action blows up as we approach degeneracy 
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For large fi, and using the approximate expression ( p2|) for the error function in (filf ) we 
have 

3 5 

which is formally the same expression as (|33|). 

In the thin wall approximation, we can estimate /i in terms of V and the bubble 
radius. For simplicity, we shall also neglect gravitational backreaction. Denoting by Si 
the tension of the bubble wall, the radius of the bubble at the time of nucleation is given 
by @ 

fl ° = (43) 

where AV = Vp — V(ar, 4>o)- Here Vf is the potential in false vacuum. For our approxi- 
mation to be valid, Rq should be larger than the thickness of the bubble wall and smaller 
than the Hubble radius in false vacuum. Under these assumptions, the Euclidean action 
is given by |TB| 

S E « ySiifr (44) 

The derivative of Se can be expressed in terms of the slow roll parameter 



and finally, from (ffTI), we have 



fx = ^R'oVe. (45) 



Taking one more derivative of InW with respect to Iny, we find 

d 2 \nW . o 7T 2 . r , / , V 



V -^(l + 4_J^) 



d(\ny) 2 

Near y = y pea k we have, setting the first derivative of W to zero and using e < 1, 

d 2 \nW 



d(\ny)< 



~ ^Vpeak' 
peak 



From this we can estimate the dispersion in the distribution of Q, which is again approx- 
imately given by (R5|). 
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VI. CONSTRAINTS FROM CMB ANISOTROPIES 



As we have shown, given a particle physics model which leads to quasi-open infla- 
tion, we can predict the probability distribution P(fl). Of course, the model also makes 
predictions for the CMB anisotropies. Comparison of all predictions with observations 
can be used to constrain the parameters of the particle physics model. 

In an open (or quasi-open) universe, CMB anisotropies which are generated during 
inflation come in three different types. The first type corresponds to scalar fluctuations 
generated during slow roll inside the bubble, and it affects wavelengths smaller than 
the curvature scale. These are called subcurvature modes. The corresponding spectrum 
of temperature fluctuations, characterized by the multipole coefficients 1(1 + l)Cj as a 
function of I, is nearly flat for I < 100. This type of fluctuations is usually believed to 
give the dominant contribution to the observed plateau in the CMB spectrum. 

The second type of anisotropies corresponds to excitations of (p generated outside 
the bubble or during the process of tunneling and expansion of the bubble into the false 
vacuum. These are accounted for by the super curvature modes discussed in section II (see 
also ]nj). For the models we have considered, the amplitude of temperature anisotropies 
caused by supercurvature modes is a factor of order Hp/WH^o) relative to the subcur- 
vature ones |30[| . However, supercurvature modes affect only the very few first multipoles, 



and hence they cannot explain the observed flat spectrum. For that reason, the constraint 
Hp 10H((p ) is usually imposed. 

Finally, there are CMB anisotropies caused by gravity waves, which can in turn be 
decomposed into the ones generated during slow-roll and the ones caused by fluctuations 
of the bubble wall itself [^[yj. Wall fluctuations give the dominant contribution for the 



few first multipoles, but their contribution decays rapidly with I. The waves generated 
during slow roll give an approximately flat spectrum, whose amplitude is much smaller 
than that of scalar modes. 

The multipole coefficients Q for the temperature anisotropies due to subcurvature 
modes are given by PU| , P2 | 



Here, we have used the notation = H(<f) ) and the slow roll parameter e given in (|40|). 
The coefficient bi is a slowly varying function of fl which can be bounded as 1 < bi < 6 
in the range .1 < fl < 1. 

Supercurvature modes induce temperature anisotropies which for the lowest multi- 
poles can be estimated as BO, 



Df c s >£±3cr^ di («F )Df (47) 



2tt \H t 
where ~ 10~ 2 . 
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Let us now consider the contribution to CMB anisotropics from tensor modes. As 
mentioned above, for the lowest multipoles this is dominated by the domain wall fluctua- 
tions [p3 j ^0l , p2l- For simplicity, we shall consider the case of a weakly gravitating domain 
wall, satisfying GSiRo -C 1, where Si is the wall tension. Also, we shall restrict attention 
to the thin wall case. Then, the anisotropies caused by the wall fluctuations are given by 



dT - -^U(n). (48) 

For the first few multipoles, and fl in the range .1 to .5, the coefficient ci(fl) is of order 
10~ 2 [for higher multipoles, q decays very fast, scaling roughly as (1 — ft) 1 ]. 

Since Hp > (8tcG/3)AV , where AV was introduced in Eq. (ff3|), we have 



HlRl > 24nG^ « 16GH 2 T ^ ] . (49) 



From (H), and using (|6]), (g7|) and (g|), we find 




(50) 

where the coefficient K\ = 4007T 2 cf dfbf /b\ Q is plotted in Fig. 4 for various values of / and 
ft. The inequality (1501) turns out to be somewhat restrictive. 
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FIG. 4. The coefficient K\ for various values of the density parameter 
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In the model given in Eq. (|l|), the parameter e = [2A(0 O )] _1 is of order 10~ 2 . From 
Fig. 4, the coefficient K\ is never smaller than 10 -4 for Q in the range .1 to .7. Hence, we 
find 

f D s V f d s V fw- w \ , , 

As discussed in Section IV, constraints from the observed fraction of clustered matter 
imply /x < 3. On the other hand, observations of CMB anisotropies require C\ < 10~ 10 . 
Hence, we conclude that this model can only satisfy all observational constraints if CMB 
anisotropies are not completely dominated by scalar subcurvature modes. 

If the observed CMB anisotropies are due to inflation, then we should have 
D SC ,D W < D s ~ 1(T 10 , and Eq.(§T]) gives /i > 100. For such values of //, the peak 
of the probability distribution is at very low values of Q, and the corresponding fraction 
of clustered matter is unacceptably small. It is therefore unlikely that the two-field po- 
tential ([!]) can give a realistic model of open inflation which will explain both a nontrivial 
value of Q and the observed spectrum of CMB fluctuations. 

This problem disappears if the observed CMB anisotropies are due to a different 
source, such as cosmic strings or other topological defects forming at the end of inflation, 
which would also be responsible for structure formation. Also, the restriction ( j50|) will 
be less severe if the observed value of Q is larger than .7, since the coefficient K\ is then 
much smaller, or in models with a smaller slow roll parameter e % 10~ 3 . 



VII. THE COSMIC AGE COINCIDENCE 

The usual objection against models with Q < 1 is that it is hard to explain why we 
happen to live at the epoch when the curvature is about to dominate. That is, why 

to ~ t c , 

where to is the present time and t c is the time of curvature domination. Observers at 
t <C t c would find Q ~ 1, while observers at t ^> t c would find Q <ti 1. It appears that one 
needs to be lucky to live at the time when Q < 1. There is another coincidence which is 
required in open-universe models and which also calls for an explanation. Observationally 
the epoch of structure formation, when giant galaxies were assembled, is at z ~ 1 — 3, or 
tc ~ to/3 — to/8. On the other hand, the interesting range of Q for open universe models 
is 0.3 < Q < 0.9, which corresponds to z c ~ 0.1 — 2, or t c ~ 0.3£o — 0.9to- We see that 
to and t c are within one order of magnitude of one another. It is not clear why these 
seemingly unrelated times should be comparable. We could have for example to t c . In 
this Section, we shall argue that the coincidence 

tc ~ t c ~ t 
may be not as surprising as it first appears. 
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Let us begin with the coindidence to ~ t c . In models we are considering here, most 
of the volume in each quasi-open bubble is occupied by regions with small values of Q, 
corresponding to small values of t c . Mathematically this is expressed by the fact that 
the "dehumanized" probability distribution dV(Q) in Eq.(14) is peaked at Q = (for 
n > 1/2). On the other hand, the "human factor" z/(fi) suppresses all values of Q for 
which t c < to, so that curvature domination interferes with structure formation. As a 
result, the peak of the full probability distribution dV(Q) = u(Q)dV(Q) is shifted to a 
value of Q corresponding to t c ~ to- Hence, we should not be surprised that t c ~ to in 
our universe. 

It remains to be explained why we live at a time to ~ to- Clearly, to could not be 
much less than tc, so we need to explain why we do not have to 3> tc- We now recall 



Dicke's observation |34| that the time to is unlikely to be much greater than tc + t+, where 
t+ ~ 10 10 yrs is the lifetime of a typical main sequence star. Noticing that t* ~ tc, we 
conclude that the expected value of to is ~ to- 

The value of t± and tc depend only on fundamental constants and on the amplitude 
of the cosmological density fluctuations. In the models we have considered in this paper, 
where Q is the only variable parameter, these timescales are fixed and one cannot address 
the question of why they are similar. 

If cooling failure (discussed at the end of Section IV and in Appendix C indeed 
represents a barrier for effective star formation, then it adds yet another timescale which 
is comparable to the other four we have encountered in this Section. This is the time t c f 
after which collapsing gas clouds of galactic mass cannot fragment and remain pressure 
supported. This timescale is also determined by fundamental constants, so the coincidence 
of this scale with tc cannot be explained within our simple model. However, it is easy to 
generalize the model so that both f2 and variable. For instance, instead of just 

one slow roll scalar field, we can consider two of them, 

2 

V{a, 0i, fa) = V t {a) + °— {g l( j>l + g 2 <f>l). 

In this case, the two slow roll fields will have different mass inside the bubble. The duration 
of inflation and the amplitude of density perturbations are determined by the point in 
the plane (0i, fa) where the fields land after tunneling. Changing to polar coordinates on 
that plane, the number of e-foldings of inflation depends basically on the radial coordinate 
R (how far we are from the bottom of the potential). On the other hand, the amplitude 
of density perturbations depends on the effective mass along the curve described by the 
inflaton, which is determined by the angular coordinate O. 

The volume factor in the probability distribution will be the same on R = const. 
surfaces, whereas the tunneling factor will choose the direction O in which the mass 
m 2 (6) oc gi cos 2 O + g 2 sin 2 O is the lowest. In our model, cr rec oc m(Q)N(R)Mp 1 , where 
N(R) ~ GR 2 is the number of e-foldings of inflation and M p is the Planck mass. Low m 
means large to, because the smaller is a rec , the longer it takes for a perturbation to go 
nonlinear. Hence, volume and tunneling factors would choose the largest possible to- On 
the other hand, t G cannot be larger than the cooling boundary t c f. Therefore, tc ~ t c / 
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could also be explained in this model. This argument can be regarded as an explanation 
for the observed amplitude of density fluctuations Q in our universe: the value Q ~ 10~ 5 
is selected by the condition to ~ i c /-0 A detailed analysis is left for further research. 



VIII. SUMMARY AND CONCLUSIONS 

We have calculated the probability distribution for the density parameter in models 
of open inflation with variable Q. This probability is basically the product of three fac- 
tors: the "tunneling" factor, which is related to the microphysics of bubble nucleation and 
subsequent expansion; the volume factor, related to the amount of slow roll inflation un- 
dergone in different regions of the universe; and the "anthropic factor" , which determines 
the number of galaxies that will develop per unit thermalized volume. It is interesting 
that the expression for the probability ([31]) depends on the underlying particle physics 
model through a single dimensionless parameter /x, defined in Eq.(|42|). 

Taking the minimum of the slow roll potential to be at = 0, the tunneling factor 
tends to suppress large initial values of 0, favouring low values of Q. However, only those 
regions for which <ft is large enough will inflate. Hence, there will be a competition between 
volume enhancement and "tunneling" suppression. 

The most interesting situation occurs when the tunneling suppression dominates 
over the volume factor. In this case, the product of both would peak at Q = 0, and the 
anthropic factor v{Vt) becomes essential in determining the probability distribution. In 
an open universe, cosmological perturbations stop growing when the universe becomes 
curvature dominated, and for low values of Q structure formation is suppressed. The 
effect of the anthropic factor is, therefore, to shift the peak of the distribution from Q = 
to a nonzero value of Q. 

As a first approximation [Tl"|,|2"0|| , we have taken u(Q) to be proportional to the 



fraction of matter that clusters on the galactic mass scale in the entire history of a certain 
region. We have found that the peak of the distribution is given by the condition 



4 Anthropic bounds on Q have been previously discussed in Ref. [|28[|. Tegmark and Rees [28 
used the inequality to < t c j to impose a lower bound on Q. To obtain an upper bound, it has 
been argued [35,^8] that for large values of Q galaxies would be too dense and frequent stellar 
encounters would disrupt planetary orbits. To estimate the rate of encounters, the relative stellar 
velocity was taken to be the virial velocity v v i r ~ 200km/ s, resulting in a bound Q > 1CP 4 . 
However, Silk Q| has pointed out that the local velocity dispersion of stars in our galaxy is an 
order of magnitude smaler than v v t r . This gives Q > 10~ 3 , which is a rather weak constraint. 
This issue does not arise in the approach outlined in the text above, since in our case large 
values of Q are suppressed by the tunneling and volume factors in the probability. 
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where the coefficient k ~ 10 _1 is defined in (p9|). For models with fi ~ 1 (which can be 
easily constructed), the probablility distribution for the density parameter V{VL) can peak 
at values of Q such that x = (1 — 0)/0 ~ 1 (See Fig. 1). The peaks are not too sharp, 
with amplitude Ay m 1/2, or Ax ~ 5, so a range of values of Q would be measured by 
typical observers. 

The analysis we presented here demonstrates that, given a particle physics model, 
the probability distribution for Q can be unambiguously calculated from first principles. 
We can also invert this approach and use our results to exclude particle physics models 
which give the peak of the distribution at unacceptably low values of Q. This gives the 
constraint \i < 3. 

An independent constraint on the model parameters can be obtained from CMB 
observations. If the observed CMB anisotropies are to be explained within the same two- 
field model of open inflation, without adding any extra fields, then we have shown in 
Section IV that the corresponding constraint (if the observed value of Q lies in the range 
.1 to .7) is \l ^ 10 6 e 2 , where e is the slow roll parameter defined in ( fjop . Combinig both 
constraints, we obtain a bound on the slow roll parameter 

e < 1(T 3 . 

This bound is somewhat restrictive. For instance, for the simple free field model ([I]), the 
slow roll parameter is of order 10~ 2 , and so this model would contradict observations. It 
is easy, however, to generalize the slow roll potential in order to make e sufficiently small. 
If one allows some other source for CMB fluctuations (e.g., topological defects), then the 
CMB constraint is much less restrictive, and simple models of the form (|I]) are still viable. 

We have advanced anthropic arguments towards explaining the "cosmic age coin- 
cidence", that is, whether it would be surprising to find that we live at the time when 
the curvature is about to dominate. We have argued that this is not unexpected. We 
have also discussed a three-field model in which the amplitude of density fluctuations Q 
becomes a random variable. We have outlined an argument explaining the observed value 
Q ~ 10~ 5 in the framework of this model. 

While this work was being completed, Hawking and Turok [[37]], have suggested the 
possibility of creation of an open universe from nothing (see also [0). The validity of 



the instantons describing this process P5[| , and also their ability to successfully reproduce 
a sufficiently homogeneous universe, is still a matter of debate and needs further inves- 
tigation. Clearly, the analysis presented in this paper can be easily adapted to this new 
framework. 
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APPENDIX A: SIZE OF THE POPULATED UNIVERSE AND CLASSICAL 

ANISOTROPIES 



As mentioned in Section II, a quasi-open universe is formed by an ensemble of 
inflating regions of very large size compared to the curvature scale. Clearly, the central 
parts of each region will inflate longer and, will have a larger density parameter than the 
peripheric regions. Hence, the fraction of clustered matter will decrease as we move away 
from the center. Here we shall estimate the size of the populated region, which, as we 
shall see, is much smaller than the size of the inflating region. 

From (|1D, (0) and flU) we have 



d In v\ peak = 2tyG(Q/j, - 3)d<j>l ~ — d In O . (53) 

This equation gives the variation of v due to the gradients in <fr as we move away from 
a typical civilization which measures the peak value of Q. [the estimate in (^) holds 
provided that /i is not too close to 1/2, say ji > .6]. 

Taking this civilization to be located at r = 0, the gradients can be decomposed in 
multipoles. For I = 0, d(po can be found from (H). For r ~ 1 (which for low Q roughly 
corresponds to the present Hubble distance) we have <iln0 o ~ 7. Combining with (|53| ) 
we find that v changes by 

fh»»^ S Jt (54) 

over the Hubble distance. 

For X <C 1, v would not change appreciably on cosmological scales. Using (H), the 
co-moving size of the populated universe can be estimated as the distance at which v 
drops by an order of magnitude, 

r v ~X-\ 

For fi > 1/2 we need </>q 3> f 2 in order to have sufficiently long inflation. Hence we find 
that the size of the populated region is larger than the curvature scale but still much 
smaller than the size of the quasi-open island, 1 C r p C l" 1 ■ 
For X > 1 we can use @ for small r to obtain 

dln0o~7^ (55) 
o 

In that case, the size of the populated universe can be estimated as 



r p ~ 



X- 1/2 < 1, (56) 



and the human factor would drop by several orders of magnitude within our Hubble radius. 
Clearly, we should not expect to lie precisely at the center of the hospitable region, but 
rather at the outskirts, and then we would observe a large anisotropy in v around us. 
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This can be confirmed by analyzing the I > supercurvature modes. The amplitude of 
/ > modes is of order [K| 



Combining with (|53|) we have 



l 1/2 fr l . (57) 



8\nv~X 1/2 r l . 



For large X, the I = 1 anisotropy in v becomes of order one at the distance r p , as expected. 

If v is proportional to the fraction of clustered matter, as we have assumed in the 
preceding section, a large drop in this quantity is already excluded by observations [4C], 
so the constraint 

must be imposed on our model. 

This constraint is relevant to the question of classical anisotropics in a quasi-open 
universe, discussed in Ref. . To an observer living at large distances from the center of 



the island r>l the universe would look anisotropic, with d<pQ ~ 70o over the curvature 
scale around that point. For X > 1 this anisotropy would be larger than the I > 
quantum fluctuations from supercurvature modes (|57|). However, as shown above, for 
X > 1 the typical observer must be at a distance r ~ X -1 < 1 from the center of the 
island, and the arguments of Ref. HT| do not apply. Hence, even though the constraint 



(|58|) coincides with the one derived in j|T| (where a single island was considered and the 
universe was not taken to be homogeneous on very large scales), its interpretation is very 
different. It does not arise from requiring that the classical CMB anisotropy should be 
smaller than the I > 1 supercurvature anisotropy but from demanding that the factor v 
determining the density of civilizations should be isotropic around us. 

On the other hand, for the simple model (|l|), one can find a much stronger constraint 
on X by combining the bounds from the observed isotropy of the CMB discussed in 
Section VI, with the bounds on the observed fraction of clustered matter. Indeed, the 
supercurvature anisotropy can be expressed as 



Df° ~ 1(TV X, (59) 

where X was defined in (|54l) . Using the constraints form the observed fraction of clustered 
matter fi < 3 |see(|36"Dl and requiring Df c ^ 1CT 10 , this results in 

X < l(ny < 10" 3 , (60) 

a much stronger constraint than d58|). Hence, the size of the populated universe should 
be at least 10 3 times larger than the curvature scale in this model. 
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APPENDIX B: EVALUATION OF k 



As mentioned at the end of Section VI, in order predict the expected values of at 
Tcmb — 2.7K in our part of the universe, we need to know a, as well as the coefficient 

3 Tcmb 
k = — -j= — — 

0\ Z<J rec -L rec 

that relates x to y. For the temperatures we take Tcmb ~ 2.7 K and T rec = 1100 Tcmb- 
The main unknown in this coefficient is a rec . 

The value of a rec is, to a very good approximation, the same in all thermalized 
regions. Hence it can be inferred from measurements of CMB anisotropies on large angular 
scales in our observable region. Since we are interested in relatively small scales, we also 
need to make some assumptions about the power spectrum of density fluctuations. We 
shall take a scale invariant cold dark matter (CDM) adiabatic spectrum. As we shall see, 
our ability to infer the precise value of a rec will be limited by the fact that the density 
parameter Qq in the observable part of our universe is not known very precisely. Hence, we 
shall leave it as a free parameter. We emphasize that fl® is the value of the density which 
is actually realized in our universe today, and whose precise value we do not know yet. 
This should not be confused with the random variable Q which appears in the probability 
distributions, and which takes different values in different regions. 

In order to determine a rec , we note that 

a rec (R) =A-\n )a (R), (61) 

where cro{R) is the present density contrast on the relevant scale R and A(Qq) is the 
factor by which linear perturbations have grown from the time of recombination until the 
present time. In an open universe, this factor is given by |23] 



L X rec 

where Xq and x rec are the values of (1 — in our observable universe at present and 

at recombination respectively. The function / is given in (p5j). 
With this, we have 

k _ 3 /(aft) 
2\/2o~o Xq 

where we have used the fact that xT = const, in the matter era. 
The present linear density contrast o"o is given by [[22] ,f20[ 



MR 9 ai) = (c 100 r) 2 ^ 1/2 (^)- (62) 

Here Cioo ~ 2997.9 is the speed of light in units of lOO/cms -1 , 5h is the dimensionless 
amplitude at horizon crossing (which can be inferred from COBE measurements), T = Q h 
is the "shape parameter", with h the present hubble rate in units of \00kms~ 1 Mpc^ 1 (we 
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ignore the effect of baryon density in this expression for T), and K contains the information 
on the power spectrum and the length scale R ga i we are considering. 
For a scale invariant spectrum, K is given by |20| 



K(R) = / q i T 2 (q)W 2 {qRhTMpc- 1 ) dq. 
Jo 

where the transfer function T in the CDM model can be approximated as 

= ln(1+ Jf 34g) [l + 3.98g + (16.1g) 2 + (5.46g) 3 + (6.71g) 4 ]" 1/4 , 
2.34g 

and the top-hat window function W in momentum space is given by 

3 

W(u) = — (sinu — u cosu). 
u 6 

In order to find numerical estimates, we shall consider |2(| R ga i = 1 — 2 h~~ l Mpc. Roughly 
speaking, this corresponds to the scale whose baryon content collapses to form a galaxy 
with a mass comparable to that of the Milky Way. Also, by requiring that CDM pre- 
dictions correctly reproduce the statistics for galaxy distribution on scales of tens of 
megaparsecs ||43|| , the shape parameter is constrained to be in the range 



T « 0.25 ±0.05 . 
For our estimates, we shall take V = .25. With this, we find 

K(l h~ l Mpc) « 0.049, K(2h~ 1 Mpc) « 0.026 . 
For the dimensionless amplitude 5h we shall use the fitting function given by Liddle 
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6 H (Q) = (4.10 + 8.83Q - 8.50fi 2 ) 1/2 x 10~ 5 (63) 
Hence, the coefficient 

« = ? f(x ) 

2v / 2(c 100 r) 2 ir 1 /2^(fi )xo V ; 

will be sensitive to our ignorance of the value of Qq in our universe, as mentioned above. 
In Fig. 1 we plot k as a function of f2 for the two chosen values of the scale R ga i- 



APPEDIX C: EFFECTS OF COOLING FAILURE 

As mentioned at the end of Section IV, fragmentation of gas clouds will only occur 
if the cooling timescale t coo \ is smaller than the timescale needed for gravitational collapse 
Tgrav Because of this, fragmentation will be suppressed after a certain critical time i*. 
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Here we shall investigate the possibility that clouds collapsing at t > i* do not 
effectively form stars even after they eventually cool. We shall see that, as a consequence, 
the peak of the distribution will be shifted to somewhat larger values of Q. 
The density of the virialized collapsing cloud p v i r is given by [|28"| , p5 



Pv 



10 2 (Gt 2 mr )-\ 



where t V i r is the time at which the collapse occurs. The virialization temperature can be 
estimated as T vir ~ m v v 2 ir ~ m p (G 3 p vir MgY^ 3 . Here m p is the proton mass, and v vir is 
the virial velocity v„ r ~ [GM g j L) 1 / 2 , where L is the size of the collapsed object. The 
later an object collapses, the colder and rarer it will be. 

The cooling rate of a gas cloud of fixed mass depends only on its density 
and temperature, but as shown above both of these quantities are determined by t vir f\ 
The timescale needed for gravitational collapse is r grav ~ t vir Therefore, the condition 
Tcooi < Tgrav gives an upper bound i* on the time at which collapse occurs. Matter that 
clusters after that time should not contribute to the anthropic factor u(Q). 

Various cooling processes such as Bremsstrahlung and line cooling in neutral Hydro- 
gen and Helium were considered in Ref. 281. For a cloud of mass M g w 1O 12 M , cooling 



turns out to be efficient f\ for 

t < U ~ 3 • 10 10 Yr. (65) 

This value of should be taken only as indicative, since the present status of the theory 
does not allow for very precise estimates. 

From the time of recombination to the time i* fluctuations will grow by the factor 
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G.(n) = ^-/(s.), 



■rec 



whereas the critical density contrast is given by 

3 

where f(x) and g(x) are given in ( p5|) and (p7|). Following the steps that lead to Eq. (|2 
we now find 



5 Actually, the fraction of baryonic matter Xf, is also relevant for cooling. Following [28 we 
shall take ~ 0.1 

6 This upper bound on t is determined by line cooling in Helium. For M g ~ 10 12 M Q there is 
also a narrow range of time near t ~ 3 • 10 11 Yr. where cooling is again efficient due to Hydrogen 
line cooling. However, the range is very narrow and we shall disregard the galaxies which may 
form during this short late period. 
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v = erfc [nxg^x*)} 



(66) 



Noting that in the matter era [23 



t = t 



x i/2( 1+x y/2 _ sinh -V/ 2 



rec 



xlic(l + x rec ) 1 / 2 — sinh 1 xVec 



and using x rec C 1 we have 



+ (^) 



2/3 1 




X 



rec 



where 




3nt. 



) 



2/3 



A, 



rec 



(67) 



2£ aft 



Therefore, the fraction of matter that clusters on a given scale before the critical time is 
basically obtained by shifting y in Eq. ( ]2"%D by the constant kA* 



Using the values flo — -5 and h = .5 for our observable universe in order to infer k (see 
Fig. 1) and t rec »s 5.6 • 10 12 (f2o^ 2 ) _1 ^ 2 s, we have kA* »s 0.2 (as in Appendix B, we have 
used T rec = 1100 Tcmb)- The peak of the modified probability distribution is plotted in 
Fig. 3 (curve c) as a function of /i, next to the original curve a where cooling failure 
is neglected. Asymptotically both curves differ only by Ay pea k « kA*/2 = 0.1. This is 
much smaller than the width of the distribution Ay ~ 0.5, so the effect is rather small. 
In Fig. 5 (a) we plot the probability distribution as a function of x = (1 — Q)/Q (at the 
temperature T C mb = 2.7 K) for three different values of /i, without taking into account 
cooling effects. For comparison, in Fig. 5 (b) we show the modified distribution when 
matter that clumps after time £* is disregarded. 

We note that even if cooling is efficient, the density of the protogalactic cloud is likely 
to affect the number and the mass distribution of stars in the resulting galaxy. Masses 
of suitable stars should be large enough to provide the necessary luminosity and small 
enough so that the stellar lifetime is sufficient to evolve intelligent life. It is conceivable 
that the number of such stars drops with the density, in which case the upper bound on 
t should be stronger than (|65|). Again, galaxy formation is not understood to the extent 
that would allow us to estimate this upper bound with accuracy. However, since we do 
not observe many giant galaxies forming at redshifts lower than z = 2, we may consider 
as a third possibility the case where matter that clumps after the time £* ~ 3 • 10 9 Yr is 
excluded from the anthropic factor. This corresponds to ft A* w 1. Even in this extreme 
case the shift in the peak \Ay pea k\ < kA*/2 w 0.5 is of the same order of magnitude as 
the width of the distribution Ay ~ 0.5 [see fl3"5|)1. The new distribution as a function of 
x is plotted in Fig. 5 (c). 

Therefore, we find that the impact of these effects on the probability distribution 
may be significant, but not dramatic, and (0) is still valid by order of magnitude. 




v = erfc(y + kA*) 



(68) 
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x = (l-Q)/Q 
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FIG. 5. The probability distribution for is sensitive to the fact that objects which collapse 
at very late times have very low density, and therefore may be unsuitable for life. Neglecting 
these "selection" effects, frame (a) shows the probability distribution for Q for various values 
of fi (The value of Q is the one measured at the temperature Tqmb = 2.7 K). In this case, 
the anthropic factor (also shown in the plot) is just proportional to the total fraction of 

matter that clusters on the galactic mass scale in the entire history of a particular region. In 
frame (b) we disregard matter which clusters after the time when helium line cooling becomes 
inefficient, so that the collapsed galactic mass objects cannot fragment into stars. Finally, as a 
more extreme case, in frame (c) we disregard matter that clumps after the time t* « 3 • 10 9 Yr, 
since we do not see many giant galaxies forming at redshifts lower than z = 2. 
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